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I.  INTRODUCTION 

The  purpose  of  this  report  is  to  present  approximations  to  the 
standard  errors  of  moment  estimators  for  the  parameters  of  the  Rician 
distribution  and  Laha's  Bessel  distribution.  The  asymptotic  properties 
of  these  estimators  will  be  investigated  and  a comparison  to  similar 
estimates  using  maximum  likelihood  techniques  will  be  performed. 

The  findings  in  this  report  are  an  extension  of  efforts  by  this 
organization  to  develop  a method  to  be  employed  in  determining  the 
frequency  distribution  of  laser-radar  return  signals  and  to  make  infer- 
ences regarding  future  signal  data  based  on  this  characterization. 


' 


I i 


II.  BACKGROUND 

In  an  earlier  study  [ 1]  supporting  research  in  the  area  of 
laser-radar  technology,  moment  estimators  for  the  Rician  distribution  and 
Laha's  Bessel  distribution  were  developed.  The  generalized  Rician 
distribution  can  be  given  by 


(1) 


is  a random  normal  vector. 


is  the  mean  vector, 


A 


> 0 


I 


is  a positive  definite  covariance  t..uci\  1 where  is  a 

positive  constant,  l is  the  identity  matrix,  and  1 ^ is  the  modified 
Bessel  function  of  zero  order. 

2 

By  applying  the  transformation  X ■ K~,  laha's  Bessel  distribution 
was  shown  to  be 


J 


l 


I 


1 


I 


f, 

K , 


r 

. 


The  mean  and  variance  of  X were  given  by 


uj  00  - A*  + 2^q 
M,  (X)  - 4i^A2  + 44  2 


(3-a) 

(3-b) 


From  Equations  (3-a)  and  (3-b)  the  following  moment  estimators 
were  found: 


A2  - +n/2M’(X)2  - M.^ 


2 . H’OO  - A 


(X) 


(4-a) 

(4-b) 


whore  M|(X)  and  M^(X)  are  the  sample  moments  about  the  origin. 

To  use  any  estimator  rationally,  its  properties  (e.g.,  its  sampling 
distribution,  expected  value,  variance)  should  be  known.  When  the  exact 
sampling  distribution  is  not  known  but  the  estimator  involves  sums  of 
sample  values  as  in  this  case,  standard  errors  provide  a means  for 
determining  the  large-sample  variance  of  the  estimators.  The  remainder 
of  this  report  will  focus  on  the  properties  of  the  estimators  given  in 
Equations  (4-a)  and  (4-b)  . 


III.  EXPECTED  VALUE  AND  VARIANCE  OF  A2  AND  £0 

The  expected  value  and  the  variance  of  the  two  estimators  can 

/\  2 ^ 

be  approximated  by  expanding  A*“  and  ^ in  a Taylor  series  about  M|(X) 
and  M^(X)  evaluated  at  ..|(X)  and  p.J(X)  wliere 

m| (X)  - E{MJ(X)| 
n£0O  - E{M^(X){ 

In  general,  if  T is  a function  of  Mj  and  then  the  Taylor  series 
of  T expanded  around  and  Uj  is 
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I 


>1^)*  (^)M..u.  <Mi  ■ “l1  * (^)  <Mi  ‘ "2' 


"iS 


* [($)_<M1 ' “p2  + ($) 
) (Ml  " m1)(M2  " m2*  + 


(M^  - n£) 


Mj)  ^ • • • • 


Taking  the  expected  value  of  both  sides  yields 


\ dMjclMj 


cov(M|Mj|) 


(7-a) 


(7-b) 


(7-c) 


The  variances  and  covariance  of  M|  and  can  be  shown  to  be 
var(Mj)  « ± [n£  - u;2] 

var(^)  " “ - M’2  j 

cov(Mj.M^)  - i ^ 

where  ^ is  the  sample  size. 

Therefore,  the  expected  value  of  T can  be  written  as 


e{t(m[,m^)}  a e|t(m;.’)|  + ^ (4  - m;2) 

+ {tvif)  ^ ^ " u^*)  * 


5 


By  considering  only  the  first-order  terms  in  Equation  (5),  the 
variance  of  T can  also  be  approximated.  Rearranging  terms  in  Equation  (5) 
yields 


T(M|,M^)  - T(uJ.p')  * (^j  (M*  - Mj)  + (^)  (M4  - mJ)  (9) 


Squaring  both  ends  and  taking  the  expected  value  yields 

2 / \ 2 


var 


[t<h;.^)]  » (5^)  var(M”  + (3^)  var(M' 


) 


+ 2 


(4)  (§■)  -»4.S> 


(10) 


Therefore,  the  variance  of  T can  be  written  as 

"*<w]’t[(t)V“i2)  +(#)2U-“  i2) 


+ 1 


(;%)  (4)  • «i“2> 


(ID 


:2 


Since  the  estimators  for  both  A and  are  functions  of  Mj  and  , 
using  Equations  (8)  and  (11)  the  expected  value  and  the  variance  of 

A ^ A 

A and  Vq  are 


E{A  } A2  + 7 ”■  7172  2m1m3  4 " 4 u2 

,Z  , \ U J 


2 1 (2,ii  -“2) 


(12-a) 


var |A“ } 


',(2Ml  - nj) 


,2 

f A *-*A  liO 

6<“2"4-l  + 4 “ 4 - 2" 


1m3 


( 12-b) 
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Using  Equation  (12)  and  substituting  the  sample  moments  for  the 
population  moments  will  yield  approximat ions  for  any  sample  taken. 
Equation  (12)  can  also  be  expressed  in  terms  of  the  population  parameters 
2 

A and  » . That  is, 


E{?[  - A2  - 


r{A2}  - 


r 8'  2 l 3 

1 4.  _2  _o_ 

.2  ■ .4  “6 

AAA 

2 1 

*2  *2  J 


zv  \|*  4\i  8s.  A * A 

E{'0l-'0  + ^ — + 


V _ „ h<L 

■ + 2 + A 

.0  A A 


+ 198 


(13-a) 


( 1 3— b) 


(13-c) 


( 1 3-d) 


As  can  be  seen  from  Equation  (13),  A*  and  are  consistent  estima- 
tors  in  that 


a)  The  expected  value  of  the  estimator  approaches  the  true  value 

as  •,  ► 

b)  The  variance  of  the  estimator  decreases  as  , 


\\ 


i 


7 


The  estimators  are  not  unbiased  for  small  sample  sizes;  however, 
the  bias  is  negligible  for  sample  sizes  greater  than  1500  tor  most  param- 
eter values  (Tables  1 and  2). 


! 


IV.  ASYMPTOTIC  SAMPLING  DISTRIBUTION  OF  A2  AND  £0 

For  the  normal  model , > ^ ( indicator  or  skewness)  is  equal  to  zero 
and  (indicator  or  peakedness)  is  equal  to  three;  thus,  the  shape  fac- 
tors of  an  unknown  model  can  be  compared  with  these  to  see  whether  a 
normal  approximation  is  reasonable.  This  approach  is  taken  in  what 
follows,  to  determine  the  asymptotic  properties  of  the  first  and  second 
sample  moments  of  Laha's  Bessel  distribution. 

Define 

h ~ N^V  V 

X^  and  X,  are  independent 

then 

/~~2  2 

R « + X,  Rician  distribution 

and 

2 2 2 

R - X^  + X,  ~ Laha's  Bessel  distribution. 

2 

R can  also  be  shown  to  be  distributed  as  a constant  times  a noncentral 
chi  square  distribution 

XJ2/'t'Q~  with  v = l and  - a^  ,Q 

X^/Yq  ~ wlth  v ■ 1 and  \2  - a^/Vp 


where  the  term  in  brackets  is  a linear  combination  of  noncentral  chi 
squares  which  is  distributed  as  a noncentral  chi  square  with 


I 

H 


8 


v - 2 2 2 

al  + a2 

\-a12+a22/'l'0  (14) 

2 2 

The  first  and  second  sample  moments  are  denoted  by  M^(R  ) and  M.,(R  ), 

respectively.  The  first  and  second  population  moments  are  denoted  by 
2 2 2 

p^(R  ) and  n (R  ) , respectively.  Now 

»1<R2)  - Vi  (k\) 

m2(r2)  - V‘2  (<2.r.) 

where  (><v2!,)  and  “2  (k\)  are  the  first  and  second  population  moments 
of  the  noncentral  chi  square. 

Substituting  sample  moments  and  cumulants  the  following  is  obtained: 
M,(R2)  - M,  (VM  <15-a) 


Mi<R2)  - *0  Ml  «x) 

m2(r2>  - <02  m2  (x;2J  . 


(15-b) 


The  objective  is  to  determine  the  theoretical  expression  for  the  first 

2 2 

four  population  moments  of  the  M^(R  ) and  M2(R  ) ; that  is, 


4 

M^R2)] 

, ^[h^r2)] 

, m3[m1(R2)] 

, m4[m1(R2)] 

wi[ 

M2(R2)] 

, u3[m2(R2)] 

• 

2 

Population  Moments  of  the  First  Sample  Moment  M^(R  ) 

“i[Vr2>]  - *‘l[*0Ml  (*£)]  ■ Vlf’l  «0] 

UjfMjtR  )J  - 


(16-a) 
( 16-b) 


I 


u3[mi(R  )J  - yQ  u3|mi  (x;J]  (16-c) 

“*M>]  - *o\[»i  «J]  (16"i) 

It  can  be  seen  that  the  problem  of  determining  the  population  moments 

2 

of  the  first  sample  moment  of  the  distribution  of  R reduces  to  finding 
the  population  moments  of  the  first  sample  moment  of  a noncentral  chi 
square  distribution. 

The  first  sample  moment  of  a noncentral  chi  square  distribution  is 

V*'2>  -*(£  XUvA) 

U=i 

and  is  distributed  as  a constant  (1/r,)  times  a noncentral  chi  square  with 

•k  * 

v = T|V  degrees  of  freedom  and  noncentrality  parameter  \ therefore, 


4 

'i»Z>] 

= 0l‘l[l 

v*,2>]  ■ 

!o  (> 

h \ 

■2  ) 

V \ / 

(17-a) 

>l(R2)] 

2 

= +0 

[VX'*>] 

♦ 2 

= “V  U2 
n 

(x;v) 

(17 -b) 

“it 

’i(r2»] 

3 

* *0  u3 

tvx'2)] 

3 

!sL 

• 3 ^3 

n 

(XA*) 

(17-c) 

A 

^0  UA 

[h,(X'2)] 

A 

vo 

" A |JA 

(x;v) 

(17-d) 

“it 

^(R2) 

1 “ ^ <V 

* * 

-f  \ ) 

“at 

^(rV 

1-^ 

* * 
(V  + 2\  ) 

m3[mi(R2)J  - — (v*  + 3K  ) 


11 


& 


4 I 4v  + 16\  . v + 4v\  + 4\ 


8 + 16A  /\|»_  4 + 8A  /tyn  + 4A  A 


- 2 |2  + 4A  /;0  + r,  + 2t|A  Aq  + 4r,A  AQ  j 

= {2)|,o4  + 4A\)3  + V + *»AV  + ^’•'o2  J 

,4[mi(r2)]  = *§  J(n  + 2)*04  + 2(n  + 2)a\03  + nA%02J.  0 

The  coefficients  of  skewness  and  kurtosis  can  now  be  determined: 


'i2[mi<*2>]  ■ 


L \4^y 


—5-  (V  + 3\) 


(V  + 2\), 


■HBBGSEB 


(2  + 2k) 


(19-a) 


72[M1(r2)]  * 32[M1<r2)] 


^4  VR  ) 


M2  M1(R> 


v0  2 
| 4(tjv  + br\k)  + (t| v + 2r\X.) 


(v  + 2k) 


( 2 2 2 2 2 
4riV  + 16n\  + ti  v + 4ri  v\  + 4n  \ 

(v  + 2k)2 


1 2 2 2 2 
8rt  + lfrj\  + 4ti  + 8r,  \ -I-  4r,  >. 

(2  + 2\)2 


3 f 2n  + 4ti\  + Ti  + 2ti  k + Ti  2k 
r.2  I (1  + k)2 


2/r>  + 4\/n  -f  1 ■+■  2k  + k 


(1  + 


ste) 


(19-b) 


As  can  be  seen  from  and  y^>  the  first  sample  moment  approaches 

a normal  distribution  as  r(  -»».  Table  3 gives  some  values  for  0 and 
0^  for  various  values  of 

An  alternative  approach  to  determining  the  population  moments  of 
2 

the  sample  moment  M^(R  ) is  by  differentiating  the  characteristic  func- 
tion of  the  mean  of  Laha's  Bessel  variate.  The  characteristic  function 
is: 


14 


TABLE  3.  VALUES  FOR 


I 


I 

. 


I 


*£(t)  - (1  - it/na)‘np  exp 


where 


(4a  (1  - it/na) 


[1] 


a 


2 % 


p-t7  p-i 


and  n ■ sample  size.  The  rth  moment  is  given  by 

drV(t)/dtr| 


^(X)  - 


t«0 


The  first  four  moments  of  M^(R2)  were  found  using  this  method  and  the 
results  agree  with  those  previously  derived. 

Population  Moments  of  the  Second  Sample  Moment  M^R2) 


' “.[*0%  (<x)]-  (<>.)] 

(20-a) 

Cv4 

PC 

'w' 

04 

■ «0J 

(20-b) 

■ *o\[M2  «0] 

(20-c) 

««[Vr2>] 

■ t0\&,  «03  • 

(20-d) 

In  order  to  evaluate  the  preceding  expressions,  the  second 
can  be  expressed  in  terms  of  k statistics: 

sample  moment 

”2  (<0 

- n : 1 k,  ) 

2 \ v,\  / 

where  the  second  sample  k statistic  is  given  by 


k2(Xv  " *7 

? ’l 


2 (c'.2(v,m)  2 "(Z  x'2t(v,\))‘ 
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Now, 


Ml[M2(R2)]  ' ♦o'  ( ' " ')  “l[k2«\)] 

(21-a) 

(21-b) 

4v*2>]  - 4A<,)} 

(21-c) 

f21-d) 

The  population  moments  of  the  1^  statistics  can  be 

expressed  in 

terms  of  the  cumulants  of  the  statistic. 

“i[vr2>]  - *o2  ( \ l)  sh  (x;!J] 

(22-a) 

“2[vr2>]-C(\ ‘)  s[k2«J] 

(22-b) 

,.3[«2(r2)]-*0<’(\  *) 

(22-c) 

*‘*[M20,i>]  ■ *08  {!LT1)  {K4[k2  (*;!>,)]  + ^ 

2 «,)])  ) 

(22-d) 

where  K is  the  population  cunmlant  of  the  k2-  statistic^. 

Expressions  for  the  first  four,  oumuiants  for  the 

statistic  can 

be  found  in  Kendall  and  Stuart  [2],  They  are 

Kl|MXi,j]  - *2  (*&) 

(23-a) 

(23-b) 

17 


32Q)  - 2)  K5(V;a)  K,(^a) 

,2(n  - n2 

8(4n2  - n + 6)  *42<x;a) 
■(2oi  - i>3 

■,(■)  - l)2 

2 


, 96(11  - 23  kj^v.o  v :.o  ,48k2  ^..o 


,('4  - 1) 


<n  - i> 


whore 


W.O 

S<0 

**<o 

S<0 

k»<\> 

v^.o 

K8(Xv,\) 


2(v  + 2\) 

8(v  ♦ 3\) 

48 (v  + 4\) 

384  (v  + 5V> 
3840  (v  + t>K) 
46080 (v  + 7\3 
645 120 (v  + 8\) 
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(23-8) 

(24-a) 
(24-b> 
(24-c) 
(24-8) 
(24-e) 
(24-f ) 
(24-g) 


V 


<25-a) 


\ - A-A0  (25-b) 

Substituting  Equations  (11)  ami  (12)  into  Equations  (10)  and  (9) 

a 1 loves  determination  of  the  moments  of  M.,(K^)  for  various  values  of  \ 

*• 

and  Figures  l and  2 show  the  coefficients  of  skewness  and  kurtosis  for 

2 2 
M.,(R  ) for  \ and  ,.  As  can  be  seen,  the  distribution  of  M^IR  3 approaches 

- 2 

normality  as  ■,  > <»  but  at  a much  slower  rate  than  N^(R*'). 

From  these  results,  it  is  concluded  that  M^(X)  and  M.,(X)  have 

as vmpotot ical  ly  normal  approximations.  Rut  M^(.X)  and  M.,(X)  asymptoti- 

callv  normal  implies  that  any  differentiable  function  of  M. (X)  and  M^(X) 

'2  ' 1 

is  also  asymptotical  ly  normal  l 3]  . Therefore,  A and  \i  , both  functions 
(differentiable)  of  Mj(X)  and  M,,(X),  are  asymptot  ical  ly  normal. 

V.  COMPARISON  OF  MOMENT  ESTIMATORS  TO  MAXIMUM  LIKELIHOOD 
ESTIMATORS 

Kst imat ion  is  a common  problem  in stat ist ics.  A meaningful est i- 
mutor  should  be  chosen  from  a class  of  estimators  having  certain  optimal 
properties  (o.g.,  unbiasness , minimum  variance,  consistency,  etc.). 

Some  of  the  better  known  estimation  techniques  are  as  follows: 

a)  Maximum  likelihood. 

b)  BAYES. 

c)  MINIMAX. 

d)  Method  vif  Moments . 

e)  Least  squares. 

Several  of  these  techniques  were  considered  for  determining  the  estimator 
2 

of  the  parameters  A and  ^ for  the  Rician  distribution  l 4],  The  method 

of  Moments  was  the  only  one  of  the  several  tried  that  admitted  to  explicit 
form.  In  previous  chapters,  the  asymptotic  properties  of  these  moments 
estimators  were  examined.  They  were  shown  to  be  asymptotically  normal 
and  consistent. 

Because  moment  estimators  are  usually  relatively  inefficient,  a 
study  was  undertaken  tv'  compare  them  tv'  the  maximum  likelihood  estimators 
•) 

(MLB)  for  A“  and  V . It  should  be  recalled  that  MLB  are  said  to  be 
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Figure  1.  Coefficients  of  skewness 


Figure  2.  Coefficients  of  kurtosi 


best  asymptotically  normal  (BAN).  The  •tati.tical  .^P^rL^uns'tere 

two  estimators  could  not  be  derived  analytically , thuslo  1 ‘ 1 

made  to  calculate  the  means  and  the  variances  of  the  estimators.  To  got 
the  MLE,  a simplex  technique  was  used  to  maximize  the  log  likelihood 
function.  The  log  likelihood  function  is, 


-» - ir  i (Ri2+A2) 


(s) 


The  moment  estimates  were  obtained  directly  using  the  equations  given 

previously.  For  each  run,  100  samples  of  size  N were  simulated.  The 

Monte  Carlo  results  are  given  in  lable  4. 

a A 

As  expected,  when  N E(A)  A,  and  E(V)  - V - VQ  for  both 

estimators.  Also  as  N > «*,  var (MLE)  < var(N0M  E) , but  not  by  much 
On  the  whole , the  agreement  between  MLE  results  and  moment  estimator 
results  is  good  and  would  imply  that  the  moment  estimators  are  almost 
as  good  as  the  MLE  in  the  large-sample  case. 

VI.  CONCLUSIONS 

This  study  has  shown  that  in  the  large-sample  case  the  moment 
estimators  are  almost  as  good  as  their  maximum  likelihood  counterparts . 
This,  along  with  the  ease  with  which  they  may  be  calculated  serves  only  to 
to  enhance  their  appeal.  Future  laser-radar  tasks  will  involve  large- 
samples  (far  in  excess  of  5000  data  points),  and  the  asymptotic  normality 

of  A2  and  V will  permit  confidence  statements  and  test  hypotheses  involv- 

2 0 

ing  A and 


TABLE  4.  MONTE  CARLO  RESULTS 


N *■  50 

N - 100 

N - 500 

A - 1.4142 
<iQ  - 0.5000 

A 

AVG  (A1  MLE 

1.4379 

1.4165 

1.4101 

VAR 

0.0184 

0.0107 

0.0015 

A 

AVG  (A)  MOM 

1.4237 

1.4135 

1.4124 

VAR 

0.0408 

0.0147 

0.0018 

A 

AVG  (V)  MLE 

0.4870 

0.4972 

0.5012 

VAR 

0.0110 

0.0067 

0.0015 

A 

AVG  (V)  MOM 

0.4945 

0.5003 

0.4977 

VAR 

0.0248 

0.0109 

0.0020 

A * 1.4142 

v0  = 1.0000 

A 

AVG  (A)  MLE 

1.4162 

1.3973 

1.4094 

VAR 

0.1129 

0.0566 

0.0079 

AVG  (A)  MOM 

1.3811 

1.3649 

1.4140 

VAR 

0.1864 

0.1215 

0.0085 

A 

AVG  (V)  MLE 

0.9825 

1.0059 

0.9990 

VAR 

0. 1025 

0.0463 

0.0111 

A 

AVG  (V)  MOM 

0.9946 

1.0167 

0.9906 

VAR 

0.1313 

0.0817 

0.0129 

A - 1.4142 
- 2.0000 

AVG  (A)  MLE 

1.1626 

1.3090 

1.3699 

VAR 

0.5271 

0.2712 

0.0727 

A 

AVG  (A)  MOM 

1.3035 

1.2566 

1.3711 

VAR 

0.5973 

0.4941 

0.1346 

A 

AVG  (V)  MLE 

2.1385 

2.0294 

2.0238 

VAR 

0.4153 

0.2592 

0.0900 

AVG  (V)  MOM 

1.9136 

1.9760 

1.9860 

VAR 

0.5623 

0.3613 

0.1239 

N - 500 


IEDE1 

N - 500 

A - 2.828A 
v|>  - 0.5000 

A 

AVC.  (A) 

2.8505 

2.8348 

2.8250 

VAR 

0.0134 

0.0064 

0.0010 

A 

AVG  (A) 

2.8483 

2.8331 

2.8254 

VAR 

0.0137 

0.0068 

0.00 10 

A 

AVG  (V) 

0.4854 

0.4931 

0.5011 

VAR 

0.0079 

0.0042 

0.0010 

> 

< 

< / 

0.4924 

0.4979 

0.4999 

VAR 

0.0103 

0.0052 

0,0012 

A - 2.8284 
\|-0  - 1.0000 

A 

AVG  (A) 

2.8421 

2.8361 

2.8229 

VAR 

0.0289 

0.0147 

0.0021 

A 

AVG  (A) 

2.8563 

2.8337 

2.8245 

VAR 

0.0320 

0.0167 

0.0025 

AVG  (V) 

0,9720 

0.9916 

1.0039 

VAR 

0.0337 

0.0183 

0.0044 

< 

< > 

0.9839 

0.9963 

0.9985 

VAR 

0.0544 

0.0266 

0.0059 

A - 2.8284 
*0  - 2.0000 

A 

AVG  (A)  MLK 

2.8716 

2.8314 

2.8192 

VAR 

0.0721 

0.0399 

0.0060 

A 

AVG  (A)  MOM 

2.8475 

2.8269 

2.8247 

VAR 

0.1432 

0.0586 

0.0073 

A 

AVG  (V)  MLK 

1.9528 

2.0133 

2.0063 

VAR 

0.1599 

0.0521 

0.0231 

A 

AVG  (V)  MOM 

1.9782 

2.0012 

1.9910 

VAR 

0.3972 

0.1737 

0.0318 
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